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The era of indirect detection searches for dark matter has begun, with the sensitiv- 
ities of gamma-ray detectors now approaching the parameter space relevant for weakly 
interacting massive particles. In particular, gamma ray lines would be smoking gun sig- 
natures of dark matter annihilation, although they are typically suppressed compared 
to the continuum. In this paper, we pay particular attention to the 1-loop continuum 
generated together with the gamma-ray lines and investigate under which conditions a 
dark matter model can naturally lead to a line signal that is relatively enhanced. We 
study generic classes of models in which DM is a fermion that annihilates through an 
s-channel mediator which is either a vector or scalar and identify the coupling and mass 
conditions under which large line signals occur. We focus on the "forbidden channel 
mechanism" advocated a few years ago in the "Higgs in space" scenario for which tree 
level annihilation is kinematically forbidden today. Detailed calculations of all 1-loop 
annihilation channels are provided. We single out very simple models with a large line 
over continuum ratio and present general predictions for a large range of WIMP masses 
that are relevant not only for Fermi and Hess II but also for the next generation of 
telescopes such as CTA and Gamma-400. Constraints from the relic abundance, direct 
detection and collider bounds are also discussed. 



1 Introduction 



If dark matter (DM) is a thermally produced Weakly Interacting Massive Particle (WIMP), 
there is a lower bound on its total annihilation cross section to avoid overabundance of 
DM in the universe: independently of the details of thermal freeze-out and whether DM is 
"symmetric" or not, its thermally averaged annihilation cross section at the time of freeze-out 
(av) should be greater than order 1 pb ~ 1CT 26 cm 3 /s. Assuming one saturates this bound 
(or in other words that the dark matter is a thermal relic), one can classify the prospects for 
the indirect detection of dark matter based on whether the annihilation is s-wave, (av) oc 1, 
or p-wave, (av) oc v 2 , where v is the relative velocity of the annihilating particles. For 
a p-wave suppressed cross section, the prospects of experimentally observing annihilation 
products are poor, since v has fallen from 0(1) at the time of freeze-out to ~ 10~ 3 in our 
galaxy today. On the other hand, an s-wave annihilator has essentially the same (av) today 
as in the early universe, and indirect detection experiments are bearing down on the correct 
range. In fact, for WIMP masses below ~ 50 GeV, the thermal freeze-out cross section is 
narrowly excluded by null results from the Fermi LAT [l}[2]. 

In this work, we explore a class of models where annihilations in the early universe and 
annihilations today are controlled by different processes and thus are naturally of different 
sizes. Such a situation is motivated in models where the DM couples most significantly to 
particles with heavy masses. For instance, in scenarios of composite Higgs and top quarks, 
the new physics sector responsible for EW symmetry breaking only has sizable couplings with 
heavy states. If the DM is part of this sector, annihilation into light Standard Model (SM) 
states is often suppressed. If the mass difference between DM and the new heavy particles 
(which we generically denote as ip) is l ess than a few tens of GeV, the DM kinetic energy in 
the early universe allows tree level annihilation into ijj pairs and the correct relic abundance 
can be obtained for couplings of order 1. On the other hand, in this "forbidden channel" 
scenario, tree level annihilation into ip pairs is kinematically forbidden in the current epocbQ 

We proposed such a scenario a few years ago in the "Higgs in Space" model [3] , which was 
inspired by models of (partial) compositeness where the top quark is the only SM particle 
with sizable couplings to the new physics sector [4j-[7]. In a simple realization we explored a 
Z' vector mediator which had large couplings to DM and to the top quark but suppressed 
couplings to the other SM particles. Thermal freeze-out in the early universe is controlled by 
DM annihilation into top pairs, and for DM masses slightly below the top mass (~ 150 GeV), 
one can obtain the correct relic abundance while having very suppressed annihilation into 
the three-body final state tWb in our galaxy. Annihilation can be dominated by one loop 
channels, and this scenario can naturally lead to large gamma ray line signatures relative to 
the continuum. 

In this work, we consider a larger class of scenarios in the same category. In all of the 
theories, we assume the DM is a Dirac fermion whose primary interactions with the SM are 
through a scalar or vector s-channel mediator, leading to annihilations of the form shown in 
Fig. [I] While there has recently been considerable interest in similar theories [8] , we flesh out 
the discussion in a number of important ways. In particular, we demonstrate that it is not 
sufficient to suppress tree level annihilations; one must also take a close look at the one-loop 

1 Another interesting aspect of this class of models is that direct detection constraints are often naturally 
eluded. 
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Figure 1: General topology studied in this paper. DM (y) annihilates through an s-channel scalar 
3> or vector resonance and loop of new fermions tp. 

contributions to the continuum gamma rays as well. We carry out a thorough investigation 
of this question by studying representative examples. 

The plan is the following: in Section [2] we illustrate the idea with a very simple model 
where DM annihilates via a scalar $ into a heavy fermion ip. We explore different choices of 
couplings and gauge representations for ip. We compute in detail all one- loop annihilation 
channels and determine the line-to- continuum ratio aimc/crcontmuum, which characterizes the 
strength of the line signal relative to the DM-induced continuum background, and provides 
a useful analysis handle on any putative line signal j9j. Section [3] discusses the vector 
mediator case, distinguishing vector and axial couplings. If DM has axial couplings, both 
vector and scalar s-channel exchanges can play an important role. We compute all one- 
loop contributions to the continuum for various assumptions on the nature of the couplings 
(scalar, pseudoscalar, vector, axial). We summarize our conclusions in Section |4j In the 
appendices, we collect the expressions describing all of the one-loop annihilation processes. 
These results go somewhat beyond the immediate needs of the current work, and are easily 
adapted to any WIMP model whose interactions with the SM are primarily through an 
s-channel mediator. 

2 Scalar Mediator Model 

As a first example illustrating the forbidden channel scenario, we consider a model with 
a minimal dark sector: in addition to the (SM singlet) Dirac fermion DM, u, we include 
another Dirac fermion ip with vector-like SM gauge interactions which acts as the portal 
to the SM. v and ip can interact by exchanging a real scalar singlet $. We will eventually 
consider various possibilities for the gauge quantum numbers of ip, but begin by assuming 
that it is is a color singlet and a doublet under SU{2) carrying the opposite hypercharge 
Y = 1/2 to the SM lepton doublet. In that case, the relevant interactions are given by 

£ D -V(y s vq> + iy^ 5 )is® - ^(y^ + iy^ 5 )^ - y H ^Hu + h.c (1) 

where H is the SM Higgs doublet and we have assumed as usual that a dark symmetry 
(parity or U(1)d) forbids v from mixing with the SM neutrinos via interactions such as 
LHv. The last term is responsible for the decay of ip — > hv (where h is the SM Higgs 
boson), with yH chosen to be small enough to avoid large ip° — v mixing (which can induce a 
dangerous Z-v-v coupling excluded by direct DM searches) but large enough for ip to decay 
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Figure 2: Representative Feynman diagrams for relevant one loop processes in the scalar 
model with an SU(2) doublet ip. The last two diagrams, corresponding to annihilation into 
7$ and ZQ, vanish due to angular momentum and CP conservation respectively. 



before BBN. For this range of values, it is safely small enough as to be roughly irrelevant for 
the remainder of our discussion. 

In addition, there are interactions such as |if| 2 $ which after electroweak breaking induce 
mixing between $ and the SM Higgs. We characterize the degree of mixing as sin a, with the 
unmixed limit obtained as sin a — > 0. A small amount of such mixing is welcome to allow $ 
to decay, and as we shall see below can help arrange for the correct v relic density, but a large 
degree of mixing would lead to a large direct detection rate which could be problematic. We 
will discuss ranges of sin a together with the relic density in more detail, below. 

The one-loop annihilation channels mediated by $ are shown in Fig. |2l including po- 
tentially three lines (corresponding to 77, 7Z, and 7$) at distinct energiesr] However, in 
practice the 7$ (and Z<3>) channels mediated by an s-channel scalar vanish because they do 
not conserve angular momentum or CP. All of the corresponding cross sections scale as 

M«((l&)V + 4( 1 &) 3 ) (2) 

which in the non-relativistic limit is p-wave suppressed for purely scalar couplings. We 
therefore concentrate on pseudo-scalar couplings of $ to DM to be able to observe a sizable 
7-ray line signal (in which case the elastic scattering via Higgs exchange after h — $ mixing is 
velocity suppressed). On the other hand, large line signals are not predicated on the nature 
of the $ — ip coupling. Detailed expressions in the v — > limit are given in the appendices. In 



2 The presence of a forest of multiple lines is a generic feature arising from SU(2) x U(l) gauge invariance 
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Table 1: Summary of the one-loop DM annihilation channels in scalar exchange models and 
typical values of the line/continuum ratio for the three charge assignments discussed in the 
text. 



addition, the final states consisting of a vector and a scalar particle vanish for the s-channel 
scalar exchange. A summary of the various channels for the scalar (S) and pseudo-scalar 
(P) $ — ip couplings is shown in Table [T] It is worth noting that the t-channel tree level 
annihilation into is p-wave suppressed in all cases, and therefore annihilation into is 
controlled by the one-loop contribution. 

As alluded to above, there is a wider variety of SM SU(2) x £7(1) gauge charges for ip 
which can result in a viable line signal. We study the particular cases of SM representations 
(£17(3), SU(2), 17(1)): 

• (1,2, 1/2): ^i/2 = ^ + ,^ ); 

• (1,2,-3/2): ^-3/2 = ^-^-); 

• (1, 1,-1): ^_i = <0- 

In terms of the dark symmetry, in the case ip = ipi/2, ip is charged under the dark symmetry 
as well and decays into vH\ when ip = ip-3/2, ip is uncharged and can decay into Pie R \ 
and when ip = ip-i, we also assume it is uncharged and decays into SM states via a tiny 
mixing with the SM right-handed lepton (note that the last term in Eq. ([I]) is only present 
for if> = ipi/2 )■ In all cases, these couplings allowing ip to decay can be chosen small enough 
that they do not play much role in dark matter searches or in determining the relic density. 

In Figures [3] and |1J we show the (complete one loop) cross sections for the cases where ip 
has scalar (dotted) or pseudo scalar (solid) coupling. We illustrate the resonant dependence 
on M$ for fixed DM mass, the dependence on M^, the line cross section to continuum cross 
section, crime/ Continuum, and contours of fixed cross sections the M v -M$> parameter plane. For 
"P — ipi/2, the WW cross section is large relative to the 77 or 7Z channels and the line 
to continuum ratio (7 Wnej (J continuum IS typically around 0.02 to 0.04 (still larger than what a 
random WIMP model with tree level annihilation would predict [9j). Much larger ratios can 
be obtained for the two other cases. For ip = ip-3/2, the ratio goes up to ~ 1, with both ip~ 
and ip contributing to the line signal. When ip is an SU(2) singlet, ip = ip-i, there is no 
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M v =135GeV, 1^=300 GeV, Y=l/2 doublet 



M^,=135GeV, M^OOGeV, Y= 1/2 doublet 





Figure 3: Gamma ray line signal and continuum strength for scalar exchange models: top: the 
(1, 2, 1/2) case; bottom: the (1, 2, -3/2) case. 
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Figure 4: Gamma ray line signal and continuum strength for scalar exchange models in the (1, 1, - 
case. In the last two plots, we added the contours corresponding to correct relic abundance. 



6 




Figure 5: Examples gamma-ray spectra for three choices of parameters leading to the correct DM 
relic abundance and for two choices of DM profiles. 



annihilation into WW, with the line to continuum ratio jumping to values as large as ~ 30. 
We illustrate the spectrum of a few promising cases in Figure |5j 

Since all of the channels experience the same resonant enhancement, the ratios do not 
depend on M$ (though as shown, together with the DM mass, it does control the absolute 
values of the cross sections). Increasing leads to the fall-off in the absolute rates shown 
on the figures. We have fixed the DM mass to 135 GeV for illustration, inspired by the 



features in the Fermi LAT data from the galactic center 11 13 ; the values of the ratios 



are not very sensitive to this choice of mass and one can easily realize strong lines at other 
energies. 

This class of models provides simple examples of potentially large gamma ray line sig- 
natures. It is also easy to obtain the correct relic abundance in this setup. For non-zero 
sin a, the h-<$> mixing produces resonant enhancement of the annihilation cross section for 
M v ~ M h /2 ~ 62 GeV, and for M v ~ Af*/2. The width of $ and its coupling to SM final 
states are both controlled by sin a. The resulting relic density is a complicated function of 
the DM mass, as illustrated in Fig. [6} Typically, there will be regions around both the h 
and $ resonances with the correct relic abundance. Obviously, the regions around the $ 
resonance will also lead to large gamma ray line features. Since the ip mass is slightly above 
the DM mass, one could also invoke coannihilation between v and ip. Similarly, when ip 
decays into SM states and it is slightly heavier than v, vv — > ipip annihilations are efficient 
in the early universe, although forbidden today. This is the forbidden channel idea proposed 
in [3], where the role of ip was played by the top quark. 

We close this section with some brief discussion of collider and precision measurements. 
Since all choices of ip are vector-like, with at most extremely tiny mixing with SM fermions, 
constraints to precision electroweak observables 14 are expected to be much weaker than 
direct constraints from colliders. Since the ip states contain charged particles which can 
decay into leptons, low-background final states at the LHC are possible; however, these are 
balanced by very small electroweak production cross sections. As a result, current experi- 
mental constraints are rather weak, and new states as light as a few hundreds of GeV are 
allowed. We expect that further high energy running of the LHC should eventually cover a 
good portion of the parameter space. 
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Figure 6: Contours of correct relic abundance in the M v -M$ plane from annihilations into ipif) and 
(where allowed) and SM states via h — $ mixing. 

3 Vector Mediator Model 

The case of an s-channel vector mediator builds on what we have learned from the scalar 
model in the previous section. The SM gauge sector is extended by a £7(1)' gauge symmetry 
which is broken by the VEV of a complex scalar, resulting in a spectrum consisting of a 
massive vector boson Z' and the associated Higgs scalar $. Their masses can be related to 
the [/(l)'-breaking VEV, 

M z , = g ziqq> ($), (3) 
M$ = y/E ($). (4) 



where qq> is the charge of the complex scalar, gz> is the U(l)' gauge coupling, and a is the 
scalar quartic, which we will fix to 1 to illustrate the typical physics. The phenomenology as- 
sociated with $ exchange is very similar to the scalar model of the previous section, including 
the possibility that it mixes with the SM Higgs h through couplings such as |/7| 2 |$| 2 . 

The DM is charged under U(l)', as are the if) fermions. As in the case of the scalar 
exchange models, there are a variety of viable choices for the representations of if) under the 
SM. In practice, we choose if) to be an SU(2) singlet carrying one unit of hypercharge. In 
principle, nothing determines the £77(3) representation of if>; in practice we choose them to 



be uncolored in order to evade large annihilation into gg 15 , 16 as well as the LHC bound 



on t'-like objects of about 450 GeV 17, 18 . We leave the case where ip is colored (and can 



mix) with the right-handed top quark for separate work 19 . 
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3.1 Mass Structure versus Couplings 

We will find that the one-loop cross sections in this setup are very sensitive to the vector 
or axial nature of the couplings, as summarized in Table [2] If all couplings in the loop are 
vector-like, there is no one-loop annihilation into two gauge bosons. In practice, an odd 
number of axial vector couplings is required in the loop for the annihilation cross section to 
survive in the non-relativistic limit. To avoid SM anomalies, we take ip to have vector-like 
couplings with respect to the SM gauge group, therefore, strong line signals require axial 
vector couplings of the Z' to v or if) (or both). This implies that both states are actually 
hybrids of underlying states with different U(l)' charges, married together by ($). 

In order to avoid mixed [/(l)'-SM anomalies, we engineer this by including vector-like 
pairs of v with charges q Ul and q V2 = q Ul — g$ and another vector-like pair of if) fields with 
charges q^ and q^ 2 = q^ — g$. In principle, the right- and left-chiral members of each pair 
could be married by a vector-like mass term, but we neglect these terms, which would induce 
dependence on a mixing angle, for simplicity 

These charge assignments insure that Yukawa interactions with $ are allowed, and after 
obtaining a VEV, result in mass terms for the two states in each pair, 

M^ = Aj($) , ilfy = A*($>, (5) 
M„ = A£<$) , M vl = \ s v , {<$>). (6) 

where ifjR = if>2R, = i>iL-, ip'n = ipiR and if)' L = ip 2 L- These masses are bounded by 
perturbativity of the couplings to be less than about < 47r($). 

This construction guarantees the presence of an axial-vector interaction. For models 
with purely vector-like interactions, we can choose q Ul = q V2 and q^ = q^ 2 , which forbids 
interactions with $, but allows us to write down a gauge invariant mass for the field of 
interest. In that case, the heavier element of the pair is superfluous. 

The most important terms in the Lagrangian read: 

C D i vi$vi + i v 2 ll)v 2 + i ipi$ipi + i if) 2 $ip 2 (7) 
- (vi(Vv* + Wu^ 5 )v2$ + h.c) - (^(4* + iy^ 5 )ip 2 $ + h.c) , 



where 



DyU = d^v - ig z <q v Z'^ v , (8) 
D^ip = d^ip - igz'q^Z' ip - igzZ^ if> - iQfeA^ if) , (9) 



for each and ipxp, resulting in equal vector couplings and opposite axial-vector couplings 
of the mass eigenstates if) and ip' to Z' as: 



£ ^ 0*'(g*x-y ^Z'^-if^Z^if) (10) 



(ii) 



5 As such terms are technically natural, there is no fine-tuning associated with this choice. 
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Figure 7: Representative one-loop annihilation channels in the vector model. The status of each 
channel depending on the nature of the couplings is summarized in Table [2j 
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Table 2: One-loop DM annihilation channels in U(l)' model due to a single fermion i/i running 
in the loop, where if) carries hypercharge 1 but is SU(2) singlet, for different combinations of 
vector (V) and axial (A) couplings. Unlisted combinations lead to vanishing cross sections in the 
non-relativistic limit. 



This module allows for tree level annihilation into Z'Z', Z'§ and ipip, all of which would 
lead to a large gamma-ray continuum. These potentially over-whelming contributions may 
be suppressed by imposing M v < M z > (i.e. A^ < gz'Q^), M v < (M$ + M z >)/2, and A^ < A5. 
Annihilation into <3>$ is p-wave suppressed, and thus there is no particular need to impose a 
relation between M v and M$. 

Representative one- loop annihilation diagrams are shown in Fig. [7| For a vector s-channel 
mediator, the 77 annihilation channel is non-zero only if DM has axial couplings (in which 
case $- mediated diagrams also contribute; see the appendices for details). This is one place 
in which our extended framework captures physics beyond what was studied in Ref. [3], 
where the DM was assumed to have vector-like couplings, and the 77 channel vanished. In 
addition, there are potentially gamma-ray line signals at various energies from annihilation 
into jZ, jZ' and/or 7$. The V<& final states proceed via s-channel vector exchange. 
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We classify different scenarios based on whether interactions of v and if) with the Z' are 
vector-like (V), axial- vector (A), or both (V+A). The various annihilation channels are listed 
along with their velocity-dependence, in Table [2j To recap the generic features: 

• One-loop annihilations mediated by an s-channel scalar require pseudoscalar couplings 
of the DM. 

• One-loop annihilations into 77 mediated by an s-channel vector are possible only for 
axial couplings of both the DM and if). 

• One-loop annihilations into require s-channel scalar exchange with pseudoscalar 
coupling of DM and scalar coupling of if>. 

• The jV channel mediated by an s-channel scalar requires vector coupling of if) to V. 

• The 7$ channel requires vector couplings of both v and if) to the Z' . 

The table displays the one-loop annihilation into Z'Z' due to the if> loops. However, while 
we have explicitly constructed a theory that is free of SM or mixed [/(l)'-SM anomalies, 
we have not been careful to impose cancellation of the U(l)' 3 anomaly. Cancellation of 
such anomalies will generically require the addition of fermions charged under U(l)' (but 
not the SM), and thus the rate into Z'Z' is somewhat more sensitive to the details of the 
UV completion. Besides, this one-loop channel, as well as Z'<3> and $$, requires additional 
box diagrams that are not considered in this work. We will always work in the kinematical 
regime where this channel is absent. 

3.2 Vector DM Couplings 

When the dark matter has vector-like coupling to the Z', a strong line signal requires that 
if) have axial-vector, or a combination of vector and axial vector couplings. In this case, 
there are three distinct lines: jZ, jZ', and 7$. We illustrate this case in Fig. [8] for different 
configurations of parameters. Typically, the jZ' and 7$ cross sections are much larger than 
jZ, with the 7$ line very large for a wide range of Z' masses. 

We show two representative gamma ray spectra in Fig. [9] obtained for M v = 300 GeV. 
In the first plot, there are two close-by lines, one at an energy of E 1 = 180 GeV from jZ' 
{a 1 z l v — 4.7 x 10~ 29 cm 3 s^ 1 ), and one from 7$ at = 223 GeV with with 07$ v = 
1.5 x 10 -28 cm 3 s~ 1 . The continuum is rather suppressed and dominated by Z$ annihilation 
with ov = 8 x 10~ 28 cm 3 s -1 . In the second plot, the mass of Z' is too large for the 7Z' line 
to be observable [a^z 1 v = 1.1 x 10~ 33 cm 3 s -1 ) but there is a striking line at E 1 = 135 GeV 
and with <r 7 $ v = 7.4 x 10~ 29 cm 3 s _1 . The continuum dominated by the Z$ channel has 
a m v = 3.9x 10~ 28 cm 3 s _1 . 

3.3 Axial- Vector DM Couplings 

When the DM has Z' axial vector couplings, it opens the door for the 77 channel and in 
addition there are contributions from $ exchange to the annihilation channels. We illustrate 
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Figure 8: One- loop annihilation cross sections and ratio to continuum for DM with vector cou- 
plings, for a variety of parameter choices as indicated. 
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Figure 9: Gamma-ray spectrum obtained from vv annihilation in the galactic center with M u = 
300 GeV and other parameters as indicated. 



this case in Figures 10, where we have chosen to switch off the pseudo-scalar coupling of 



the DM, and thus all diagrams mediated by s-channel scalar exchange vanish. 

In Figure 10 , we show the dependence of the cross sections on Mz>, M„, and M$. Because 
all three of these parameters are proportional to (<&), we vary them in a controlled way, 
exploring successively the variation with changing ($) and all couplings fixed, variation of 
My obtained by varying with all other parameters fixed, and variation of M$ obtained 
by varying with all other parameters fixed. From the figures, we see that either the 77 
or 7$ may end up as the strongest of the line signals (with 77 vanishing when the Z' goes 
on-shell, as required by the Landau- Yang theorem [20]), whereas the one loop continuum is 
rather anemic and is almost entirely made up of the Z$ channel. 

We show a representative two-peak gamma ray spectrum in Figure 11 for a case with 
M v = 330 GeV. For these parameters, the ^yZ' line is somewhat too faint to be observed 
compared to the Z$> continuum {a^z> v = lx 10 -29 cm 3 s -1 and <Jz<$> v = 4.5 x 10 -27 cm 3 s _1 ), 
but there are two large lines from the 7$ and 77 channels at energies of E 1 = 155 and 330 

GeV, respectively, with cross sections <r 7 $ v = 8.5 x 10~ 28 cm 3 s _1 and <r 77 v = 1.2 X 10~ 28 

3 — 1 
cm s . 



3.4 A Purely Vector-Like Model 

There is one further ingredient that is worthwhile to consider adding on top of the minimal 
framework we have introduced so far. In theories where the v and tp fields are composites, 
the most natural situation is one with only vector-like couplings. As shown in Table |2j 
this is not expected to lead to any observable line signals. However, such models often also 
predict additional scalars which are unrelated to the £7(1)' breaking, and are singlets under 
the SM, but have scalar interactions with the ip fields. For example, composite Higgs models 



based on the SO(6) — > SO (5) symmetry-breaking contain such a state S 22,21 . Such an 
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8VvZ— SAvZ'- l/2=g V y Z >=&; i j / >2!=*gj al/ z— g Al p'Z' 
M v =0.75 M z ., =0.85 M z ., M^.=0.8 My, M^My 



SvvZ'- SavZ'- ^-gv^z'-Sviii'z— 8At1!Z— 8Ai^'Z' 
M z .= 1000 GeV, M, = 0.85 M z ., M T ,=0.8 M z ., M t =0.5 M z , 




400 600 800 1000 1200 1400 
M z , [GeV] 
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M v [GeV] 



Svv Z ,_ S Av z'~ 1_ 8v^ Z ' _ Svhi Z ' _ 8 A ti)Z' _ Satp'Z' 
1^=330 GeV, M t =480 GeV, M z , = 600 GeV, M^=l.l M . 




300 350 400 450 



Figure 10: Cross sections for models where the DM has axial vector couplings, as a function of: a) 
(<3?) (with all couplings held fixed); b) y£$/M„ (with all other parameters held fixed); c) y^/M^ 
(with all other parameters held fixed). 




Figure 11: Two-peak Gamma-ray spectrum obtained from vv annihilation in the galactic center 
with M v = 330 GeV and other parameters as shown. 
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additional state opens up the possibility of strong gamma ray lines from the 'yS channel 
mediate by Z 1 exchange. In addition, one would also expect annihilation into ZS and Z'S. 
In order to simplify the discussion, we assume that the DM coupling to S is scalar, and 
thus any additional mediation via S itself vanishes in the non-relativistic limit. This case 
is summarized in Table [3] and illustrated with some representative quantitative results in 



Figure 12 
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Table 3: One-loop DM annihilation channels in the purely vector-like U(l)' model with an addi- 
tional gauge singlet scalar particle S. 



SvvZ'= 2 . Sv^z'=2,y ,ps =1 
M v =135GeV, ]VL=600 GeV, M s =150 GeV 




200 300 400 500 600 700 800 
M z , [GeV] 



Figure 12: One-loop annihilation cross sections into "yS and ZS, for the purely vector-like Z' 
model. 
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3.5 Summary 



To summarize this section, we have seen that in a simple model with vector mediator, there 
is a vast range of possibilities for large line signals depending on the nature of the couplings 
and mass relations. A key-result is that axial couplings are required to have non-vanishing 
1-loop cross sections into two gauge bosons. Gamma ray spectra with several lines are typical 
in the presence of axial couplings. We have presented a minimal anomaly-free construction 
that realizes non-vanishing axial couplings. It relies on two heavy fermions ip and ip' with 
equal vector couplings and opposite axial couplings to Z' . 

If My = M^i, the annihilation channels into two gauge bosons, namely into 7Z, 7Z 7 , ZZ 
and ZZ 1 ', vanish as the two diagrams with respectively tp and tp' in the loop cancel exactly. 
On the other hand, for 7$ and Z$ channels, no such cancellation occurs as these channels 
do not require axial couplings. In general, and M^i are different, however there is some 
partial remnant cancellation which reduces the jZ and ZZ channels relatively to the others. 
The jZ' and ZZ' channels are not as suppressed (when kinematically allowed) because of 
the larger couplings involved, g^ t z r ~ 10gy,,z, and two possible axial vertices. 

We therefore conclude that quite generally, the Z$ and 7$ channels are the dominant 
1-loop annihilation processes. 



4 Conclusion 

Very large gamma ray line signals would motivate models where DM annihilation rates are 
enhanced by resonant effects and where the continuum is depressed because of annihilation 
into heavy channels. In this work we have explored the two main categories which fit easily 
into the framework of a renormalizable quantum field theory: vector-resonance and scalar- 
resonance models. We have introduced simple effective field theory sketches of such theories 
where the dark matter is a Dirac fermion, and performed the first complete analysis which 
includes the one loop contributions to the continuum annihilation modes. 

Both cases, with appropriate choices of couplings, can lead to large line signals as well as 
large line/continuum ratios. For s-channel vector mediators, this requires chiral couplings to 
the fermions running in the loop, whereas for scalar mediators, pseudo-scalar couplings with 
DM are necessary. The various cases are summarized in tables [T] and [2] We have restricted 
our models to those in which the fermions running in the loops are disconnected from the 
SM fermions. We explore the case in which the loop fermions mix with the top quark in a 
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Gamma ray lines are a fascinating, powerful probe of dark matter annihilations. Mapping 
out the kinds of theories which naturally produce such features is an important step in 
preparing for a discovery. 
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A Effective Vertices and Amplitudes-squared 

In the following appendices, we summarize the one-loop expressions for the effective vertices 
and amplitudes-squared needed for the calculations performed here. The generic topology 



for the loop diagrams considered here is shown in Fig. 13 We express all amplitudes in 



terms of two-point (B ) and three-point (C ) scalar integrals where: 



C = C {MlMlAMl-mlmlml), (12) 

5 (23) = B^Mlmlml), (13) 

S (13) = B {AM^mlml), (14) 

B„(12) = B^M^mlml). (15) 



In the following expressions, for the sake of simplicity, we have set all loop masses equal 
(mi = m 2 = m 3 = m/). However, to derive our results, we have computed all of the 
expressions in these appendices in terms of the general masses as depicted in Fig. [13} 
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B One-loop annihilation into 77 (and gg) 

The annihilation channel vv — > 77 can proceed via s-channel exchange of both a Z' as well 
as the scalar <3>. The effective vertex for — y ^{pijYiPi) takes the form: 

V^ 77 = Ag^ + + iCe^^ , (16) 

where the coefficients A, B and C are: 

^1,2)^(2,3)^(1,3)™/ i 2C { m } ~ M l) + !) 



2tt 2 



(17) 



R ^(1,^(2,3)^(1,3)™/ (2^0 K - M D + 1) 

B = ^Ml (18) 

r ^(1,3)^(1,2)^2,3)™/^ nq v 

C " 2^ (19) 

The couplings ^ and 2/*^ are the scalar and pseudoscalar couplings, respectively, of $ to 
the fermions running in the loop, while (o^-)) are the vector (axial-vector) couplings 
between fermion pair and gauge boson V. 

The effective vertex for Z' — y 77 can be written as: 

Vz^ 77 = , (20) 

where we have used transversality of the photons (i.e., e(pi) • p\ = tip?) • P2 — 0) as well as 
the fact that the photons are identical particles to eliminate other possible tensor structures. 
The loop coefficient D is given by: 

a z '- 

D = -JL [I + 2m}C ] . (21) 

Averaging over initial spins, summing over final spins and including a factor of 1/2 for 
identical particles in the final state, the total matrix-element-squared for v + v — y 77 is given 
by: 

\M vP ^*\ 2 = \M^\ 2 + \M Z >\ 2 + 2Re |M* • M*,| , (22) 
where the matrix-element-squared for the ^-mediated amplitude is: 

]T|.M|| = " [\A\ 2 + AMt\C\ 2 } . (23) 

The denominator factor is defined as: 

= 4M^ - Ml + iM^Ts, . (24) 
The amplitude-squared for the Z'-mediated channel is: 

= = 64 (a£ )^(4Ag-M|-) 2 

^ M|, |E Z ,| 2 
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where: 

H z ,= AMI ~ M h + iM z ,T z , . (26) 

It is interesting to note that, in the limit where the Z' goes on-shell, we replace 4M^ — > M|, 
and the amplitude-squared vanishes as required by the Landau- Yang theorem. The cross 
term between the $- and Z'-mediated processes is given by: 



2Re\M tS> - M%,\ = -Re 



Uyl vD a z v M (4M, 2 - Mf,) 



(CD* + C*D) 



(27) 



Finally, we note that we can obtain the vv — > gg channel with the simple rescaling: 

2 



'99 



a 



1 



Ots Qt 



(28) 



where Q t is the charge of the top quark in units of e (i.e., Q t — 4/3). 

C One-loop annihilation into 

The annihilation channel vv — > can proceed via both s-channel $ and Z' exchange. The 
effective vertex for the $$$ vertex takes the form: 



V. 



<j>_^(j><j> 



= A 



(29) 



where: 



A = 



im,f 
2^ 



yfjf (vtff) ( fl o(13) + 2B (23) - (2Ml + M|)C 

{ytfif ( B °( 13 ) + 25 o( 23 ) - ( 2M * - 4m / + M ^ c ° 



while the effective vertex is given by: 

The loop coefficients (B and C) in this case are: 

• Z' <f> <S> 2 /~t 

m ff-y P jfys,ff m f c ° 



B = C = 



The amplitude-squared averaged/summed over initial/final spins takes the form: 
\M uP ^\ 2 = \M$\ 2 + \M Z >\ 2 + 2Re |M* • M* z ,\ , 



(30) 



(31) 



(32) 



(33) 
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where the individual parts are: 



2Re |Af* • M* z ,\ 



fa%) MS iAi 2 

| Zj<J> I 



(34) 



MJ 



(M^(a^) 2 (5 + + C*) 

+ {v^y (B - C)(B* - C)\ - (v z v '^Ml(B - C)(B* - C* 
+16 [a z ^ M\B + C)(B* + C*) 

-8 (a%S 2 M 6 M|, (5 + C) (fi* + C*) , (35) 



2y% vV v%Ml (4M, 



m 2 : 



A(B* + C*) 



(36) 



Note that a consistent 1-loop calculation of the vv — y $$ channel requires the addition of 
box diagrams. This is beyond the scope of this paper which is focussed on the topology 



shown in Fig. 13. In all our illustrative plots, we assume M v < Ms- 



D One-loop annihilation into V7 

The annihilation vv — y V^j (where V = Z or Z') can proceed via both $ and Z' s-channel 
diagrams. The effective vertices take the forms: 



and: 



V 7.1 



Ot[lV 



£) e Pl^a _|_ ^ e p2ixva _|_ jp^a^pwwv _|_ Qp^PiP 2 ^ -|_ Hp y l 6 PlP2 ° ltl + Ip tJ 'e PlP2011 ' 



The loop coefficients are: 

7 v $ 

4 = ffffVsjf f ^ 1 ^2 1 _ 2j g (i3) + 2ff (23) + C Mft + 4C m^ + 2 
4tt 2 I 4M2 - 



+16C M, 4 - 8Mj C + 2ml + 1 



5 



v lf v H<ff m f ; 
2vr 2 ( M 2 - 4M^ 



My -2S (13) + 25 (23) + C M 2 + 4C m^ + 2 



-16C M, 4 + 8M 2 ( C ( My + 2m 2 f + 1 



(37) 



(38) 



(39) 



(40) 
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ytjf v lf v H- m f c ° U1 , 

2tt 2 ' 1 ' 

v] f (2B Q {l2)M 2 v {a] f z ' f + afrf f ) ( 4M * + M v) Wf v ff + 

8^ ^ M 2 _ 4M, 2 + 



M 2 - AM;, 



2B (23)M 2 UmI + ) (tffff + a // u //) 

+ 2(B 1 (12) + 1) (a faff + afjrffj) 



Ml - AMI 



A „ 2< z , v v z> , SC aJ-vf-M 2 M 2 8C a z -vJ-M 2 M 2 
+AC om 2 (a z f -vJ f - - a rrf f ) } » ^ ^ ^ ) , 



v ]f( a] ff v ff + a ff v ¥f) . 

X 



8vr 2 ^4M 2 - Ml 

[Ml (s (13) - 5 (23) - 2ComJ + 1^ + M 2 (%C Q m) - A^j , (43) 

// V // // f±JJl I £ ( 13 ) / 16M 4 + 16M 2 M 2 + M 4 \ 

4tT 2 ^4M 2 -M 2 ^) ^ ^ ' 

-2 (m£ ^4M„ 2 ^ (12) + 25 (23) + C Ml - 2C m 2 f - 
+M 2 ^-5 (12) + 5 (23) + 2C m 2 f + lj - 16C m£J 

-B 1 (12)(^Ml-4M^j (44) 



v 1 -(a v -v z '- + a z '-v v -) 



4tt 2 



(V 2 - 4M 2 ) 

(b (13) - B (23) - 2C mJ - 1^ + M 2 (sC m 2 f + ±J ^ , (45) 
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ff y ff ff ff ff' I Q i ii i ii 
4tt 2 I 4M 2 - Mi 



H = — 11 '/ 31 - L ^\ -8M 2 M V [ B (12) -350(13) + 3B (23) + 2^(12) 



^4M 2 - M 2 ) 

+C M V + 2C mJ + + M^(25 (12) - B (13) + 2^ (12)) 

+16M^S (13) + 2^(12) + C Q M v + 2C Q m) + l^j , (46) 



v J -(a v -v z '- + a z '-v v -) 
j = ' fp '/J ff 'ff if x 



4tt 2 (My - 4M 2 

My (b q {U) - 5 (23) - 2C m 2 f - lj + M 2 (^8C m 2 f + 4^ (47) 
The amplitude-squared averaged/summed over initial/final spins takes the form: 

\ M ^v-y\ 2 = \M*\ 2 + \M Z >\ 2 + 2Re |M 4 • M*,| , (48) 
where the individual parts are: 

\M,\ 2 = A{V * v /. 2 Ml U \A\ 2 + \C\ 2 (M v - AMI) 2 ) , (49) 
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(v%y \2{DE* + D*E) (4M 2 ~ My) ( 12 M 2 + M v) 



2\ 3 



+ (F>F* + D*F) [AMI - M^) 3 + (DG* + F>*G) (My - AMI) 
A(DH* + D*H)Ml (AMI ~ M vf \ D ? (4M 2 + M v) ( M v ~ 4M 2 



2^ 2 



My 



+ 



My 



2\ 2 



-(FF* + E*F) {AMI + My) (My - 4M 2 ) 



-(FG* + F*G) (4M 2 + My) (My - AMI) 2 + 8(F#* + F*#)M 2 (My - 4M 2 ) 



r2\2 



1 



+ |F| 2 (16M* + 24M 2 My + My) + -(FG* + F*G) (4M 2 - My) (M y - 4M 2 ) 

2 



+ ^|F| 2 (M 2 - 4M 2 ) 4 + i|G| 2 (My - 4M 2 ) 4 + 



4 |FfM 2 (M 2 -4M 2 ) 4 * 



+ 



/ z ,\ 2 / 2(F>F* + D*E) [AMI - M$f (Mf, - 4M 2 ) 2 

W v mi, 

{DF* + F>*F) (4M 2 - My) 3 (M|, - 4M 2 ) 2 



M 2 



M|, 






(F>G* + F>*G) (4M 2 + My) (My - 


• 4M 2 ) 2 (Mf , 


- 4M 2 ) 2 


Ml, 






|F| 2 (My - 4M 2 ) 2 (Mf, - 4M 2 ) 2 






Mf, 






(FF* + F*F) (4M 2 - My) (M 2 - 


4M 2 ) 2 (Mf, 


- 4M 2 ) 2 


Ml, 






(FG* + F*G) (4M 2 + My) (My - 


4M 2 ) 2 (Mf, 


- 4M 2 ) 2 


Ml, 






|F| 2 (My - 4M 2 ) 2 (Mf, - 4M 2 ) 2 






' Ml, 






| (FG* + F*G) (4M 2 - M 2 ) 3 (4M 2 + M 2 ) (Mf, 


- 4M 2 ) 2 



|F| 2 (My - 4M 2 ) 4 (Mf, - 4M 2 ) 2 



4Mf, 
r2^2 



|G| 2 (4M 2 + My) z [Ml - AMI) 2 (Mf, - 4M 2 



2Re\M^-M* z ,\ = Re 



AM\, 

^2/ P %M (My - 4M 2 ) 2 (4M 2 - Mf,) 



x 



Mf ,S<j>S^, 

G (4M 2 (F* + G*) + My(G* - F*) + 2F>* - 2F*) 



(50) 
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E One-loop annihilation into VV 

The annihilation channel vv — > ViV 2 (where possible final states include ZZ, Z'Z', ZZ' and 
gets contributions from both s-channel $ and Z 1 exchange. The amplitude from 
the ^-exchange diagram takes the form: 

V^viva = + Bf 2V \ + iCF*** , (51) 

while the effective vertex for the Z'-exchange diagram is given by: 

V^vj = De Pl ^ a + Ee P2 ^ a + Fp«e PlP2 ^ + Gp%e PlP2 ^ + Hp u 1 e PlP2a ^ + Ip%e PlP2av . (52) 

The loop coefficients for the ^-exchange diagram are (for simplicity, we assume V\ = Vq): 

A = ~ gJiMpM*) ( ( a //) 2 ( 3Eo(12) + 2Eo(13) + 3Eo(23) ~ 2CoM v 
+8C m 2 f - 4^ - 4M 2 (5 (12) + 5 (23) + C M 2 + 2C m 2 f - l) + 8C M 4 

+ (vj f y ^— My (£ (12) - 2fi (13) + 5 (23) + 2C M y + 8C mJ + 4) 

-8C M* + 4M 2 (3C M 2 + 2C m 2 + l)^ , (53) 

5 = ^.// m/ ( C y \ 2 t 2M 2 _ M 2\ f M 2 [ 4S ( 13 ) _ 45 (23) - 2C M y 

16tt 2 M 2 (M 2 - M y ) 2 ^ '\ V 

+4C m* + 2^ + M y (-S (13) + 5 (23) + C M 2 - 4C m 2 f - 2) + 4C M 4 

+ (^-) 2 ^—2M 2 M 2 (-S (13) + 5 (23) + 2C M y + 6C mJ + 3) 
+M 4 (B (13) - B (23) - C M y + 4ComJ + 2) - 8C M 6 

+4M 4 (2C (2M y + m J) + l) ^ , (54) 



,2' 



C = - 2^(M-X)(M + Mv-) l B ° (13 n a //) -^>K 



V 



2 



+C M 2 [a] f - v v ff ) [a] j + v v ff ) + C (vty * , (55) 
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while the loop coefficients for the Z'-exchange diagram are 
1 / ,// \ 2 / \ 2 



D = ^ (4 ((«//) + (*//) ) + 2 ^/>?/J - 2a / >?^ (13)(a^-) 

-S (23) (a^-) 2 + C M 2 (tf f - - v v fI )(a v f j + + C (u^)* , (56) 

2af '-m 2 , / / \ 2 / \ 2 

* = 4tt 2 ( Ml — My) p (13) K) - 5 °( 23 )(°//) +^(ay / -< / )(a^- + ^ / ) 

+C (^-) 2 M*) - [M v - M V ){M V + M V ) (af f {[a]-^ + (vj f )^ + 2a / > / > / ^7) 
F = 0, (58) 
G = 0, (59) 



( a &((</) 2+ (</) 2 ) +2a 
# = ^ - - 7 -2 / -(2M 2 M 2 [B (13)-B (23) + 2C M' V 



32vr 2 M 2 (Af2 - M v f 
-6C m 2 f - 3^ + M v (B (13) - B (23) - C M V + 4C m 2 f + 2) 

+M*(8C m} + 4)), (60) 



2.U;.Uf B (13) - 5 (23) + 2C M^ 



32tt 2 M 2 (M 2 - Af2 f 
-6C m 2 f -3^j+M v (B (13) - B (23) - C M V + 4C m 2 f + 2) 

+M*(8C m 2 f + 4^j . (61) 

The contributions to the total amplitude-squared are: 

8 (v® ) 2 M 2 / / 
\M*\ 2 = -^^——[8M 2 (M-M V )(M + M V )(-M V ((AB* + A*B) + 2\B\ 2 M 2 ) 
M v \h$\ z \ \ 

+2(AB* + A*B)M 2 + 2\B\ 2 M 4 + \C\ 2 M v ^j 

+ \A\ 2 (4M 4 - 4M 2 M 2 + 3M v )^j , (62) 
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^) 2 -2(^) 2 ) (2(DE* + D*E) 



-\D\ 2 - \E\ 2 ^j + 8 (v%y Ml{{DH* + D*H) - (DP + D*I) 
-{EH* + E*H) + (EI* + E*I)) + 16 (w^) 2 M^tf] 2 + |/| 2 )^ 
+M 2 U M*(-2(DE* + ((af:) 2 - 5 (v^ 

+8 (v^ 2 M 2 {-3{DH* + + (£>/* + D*I) + (£#* + £*#) 

-3(EF + £*/)) + (|D| 2 + \E\ 2 ) ([atf - 3 (y*)' 

-32 (^) V 4 (|iJ| 2 + |/| 2 )) +A( y vf D ) 2 Mt<^Ml[{DH* + D*H) 
+(EI* + + M 2 (\H\ 2 + |/| 2 )^ + |L>| 2 + \E\ 2 ^j 

+16 (afpj 2 M®My (M v - M V ){M U + M V )(-2(DE* + D*E) + |Z?| 2 + \E\ 2 ) 
-8 (af;) 2 M*M V M 2 ,(M U - M V )(M V + M V )(-2(DE* + 
+ |D| 2 + |£| 2 )) (63) 



2i?e |M* • M* z , 



Re 



32a z ' p y* uP M*C(M 2 - M v ) (4M 2 - M 2 ,) (D* - E*) 
M 2 ,T<s> (IV)* 



(64) 



F One-loop annihilation into 7$ 



The annihilation mode vv — > 7$ can only proceed via s-channel Z' exchange. The effective 
vertex for the Z'7<3> interaction takes the form: 



Ag a>1 + iBe PlP2a » , 



(65) 



where the loop coefficients are 
A = - 



v Jf<ff v ff m f 



4tt 2 ( 4M2 - M 2 



8M 2 S (13) - 5 (23) - C M£ + 2C m z f + 1 



+16C M 4 + M| (c (m| - Am) 



- 2 



B = 



7 <E> Z' 
V .ff y p,ff V ff 

2tt 2 



m/C 



(66) 
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The amplitude-squared is: 



\M\ = 



2 V, 



2_ 
vv 



4|S Z ,| 2 



4 |A| 2 + (4M 2 - M|) 2 |5| 2 (67) 



G One-loop annihilation into V$> 

As for the 7$ channel, the annihilation into V<& (where V — Z, Z') can only proceed via an 
s-channel Z'. The effective vertex is given by: 

V%+v = A 9 a " + Bfirf + Cfirt + iDe^ a . (68) 
The loop coefficients for this process are: 

A — J J 1 x 



I . ( -2ia [ \M* ■■- M* ) + { Mr - H/2 ) +Mi 



y* Jf af f (-4MZ (-2M V ^B (12) + B (13) + 25 (23) - C M| + 4C m} - 2^j 
+M| (-AB (12) - 2B (13) - 2£? (23) + C M| - 8C m 2 f + 4^ + C M^j 
+ (M V - M $ )(M y + M*) (m| ^2B (12) + 2S (13) - C M| + 4C mJ - 2^ 
-2M 2 ^0(13) + 5 (23) + 2C m) - l^j + C M^j - 16M^2S (12) + 25 (23) 
+C (^My + Ml + 4mJ^ - 2^ + 64C M^ + yf J} vf f (-AM^ (W 2 (b (12) 
+B (13) - 25 (23) - C M| + 4C mJ + 2^ + M% ^25 (13) - 25 (23) - 3C M| 
+8C m 2 f + 4^ + C M^ + (M y - M*)(M V + M*) (^2M 2 ^S (12) - 5 (23) 
-C M| + 2C m 2 + 1^ + C M^ + M| (c M| - 4C mJ - 2^ 
+1 6M* ^2S (13) - 25 (23) - C (m£ + 3M| - 4mJ^ + 2^ + 64C M^ ^ , 



(69) 
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R tff M v m M,fi a ff + yl s f z ff) 

£> = 2 2 

-16M* ^2B (12) - 4B„(13) + 25 (23) + C ^M v + M| - 4mJ^ - 2^ 



-4M^ ( 2M F ( S (12) + S (13) - 25 (23) - 3C M| + 4C mJ + 2 J 



+Ml ^-4S (12) + 2S (13) + 2B (23) + C M| + 8C m} + 4 J + C M* 
+ (M y - M$)(M y + M*) ^M^4S (12) - 2S (13) - 2B (23) + 4C mJ + 2^ 
-Ml (-2B (12) - 250(13) + 45 (23) + C M| + 4C mJ + 2^j + C M y ^ 
+64C M^ , (70) 
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C = V J^l 

2n 2 (2M u - M v + M^f(2M v + M v + M $ ) 2 ( M 2 - (M* - 2M) 2 



ytjf a ff{ AM » ~ M v~ M l) (^-16M^2B (12) - 4S (13) + 2B (23) 
+C ^My + M| - 4mJ^ - 2^ - 4M, 2 ^2M£ ^ (12) + S (13) - 2£? (23) 
-3C M| + 4C mJ + 2^ + M| ^-45 (12) + 2S (13) + 25 (23) + C M| + 8C mJ + 4^ 
+C M v ^j + (M v - M^)(M V + M*) ^My ^45 (12) - 2S (13) - 25 (23) + 4C mJ + 2^ 
-M| ^-2S (12) - 250(13) + 45 (23) + C M 2 + 4C mJ + 2^ + C M v ^j + 64C M^j 
+ yf Jf vf f ^-32M^M v (^B (12)-3B (23)-5C Ml + 6C m} + 3^ + M|(^2£? (13) 
-25 (23) - 3C M| + 6C m 2 f + 3^ + 8M* (V* (Vb (12) - 3^(13) - 35 (23) 
-2C Ml + 6C m 2 f + 3^j+ M V M% ^6B (13) - 65 (23) - C M| + AC Q m) + 2^ 
+M| (s (13) - 5 (23) - 2C M| + 6C mJ + 3^ + C M^j 

-(M v - M^)(M V + M*) ^6^0(12) - 4B (13) - 25o(23) + 3C M| + 4C mJ + 2 
+M v Ml ^65o(12) + 4S (13) - 10£? (23) - 5C M^j + C M y + M| (c M| - 4C mJ 
-2^ + 128M^S (13) - B (23) - C (m 2 + 2M 2 - 2m£j + lj + 256C M^j ,(71) 

aXfirif 

D = ^— g x (72) 

2vr 2 ^-2M| ^4M 2 + M 2 ^ + (V 2 - 4M 2 ^ + 

^1^*^(50(12) + 5 (13) - 25 (23)) - C yf ff -vf f -(^M 2 u + M^j 

+ (aM 2 -M 2 ^ (-2B (12)yl ff af f + 2B (13)y* ff af f 

+C (AMI - *4) Kfrf'f + Vtfffj)) + CoMi(yt ff -vf f - y^-af^ (73) 
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Finally, the matrix-element-squared is given by: 



+ 



+ 



+ 



+ 



+ 



1 


'{AB* + A*B) 




VoM 2 M% 



Ml - AM 2 



-2M,f,( I.1A- h Mf ) 4 ( .Uf - n/ J ! + M. 

2 



X 



4M 2 — My + Mg ) + 



M|, ( -4M 2 — My + M, 



(AC* + A*C) 

wmJmJ, 



-2M U + My - M*)(2M„ + My - M*)(-2M„ + My + M$) x 
(2M„ + My + M*) (^4M 2 + M 2 - M 2 ^ ^(a^ (Vf, - 4M 2 



+ u 



Ml,jl6MyMl, 



8M 2 M|, 



Mf , - 4M 2 



+ u 



L 1 ' /, | I M : Uf. )-(;/; I / ) .U, 

2 



16M 4 + 8M 2 ( 5M V - Mi ) + ( M, 2 7 - M 



(5C* + 5*C) 



32M 2 M|, 



2 



-2M| ( 4M 2 + M 2 ) + ( M 2 - 4M 2 ) + Mi \ x 



.z' 



M|, - AM 2 



)(■ 



16M 4 - ( My - M, 



M 4 . ( -2Mi ( 4M 2 + My ) + ( Mv - 4M 2 ) + M, 



\B\ 



32M 2 M|, 



(«S) 2 
\C\ 2 



32M 2 Mf, 



-2M| ( 4M 2 + M 2 ) + ( M 2 - 4M 2 ) + Mi ] x 



M 2 ., - 4M 2 



2 / \ 2 

l.U 2 - M 2 + M| 



+ 



L v 



M% I -2.1/ 2 ( l.U 2 + M 2 ) + ( M 2 - 4M 2 ) + M, 



2M 2 , ( l.U 2 + My ) + [ M\ 



I U ) +M|) x 



Ml, - 4M 2 



Ml -2Mi 4M 



2 / \ 2 

4M 2 + Ml - Ml 



2 t v Z>\ 



) 



M< 



-2M| I 4M 



My 



2 ' My 



M 2 - 4M 2 



My - l.U 2 



+ M| 



Mi 



(74) 
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